Computer apparatus, program and method for determining the equivalence of two algebraic functions

ABSTRACT

In a compiler or program proving or verification software code, it is often useful to analyse a computer program and one way in which this might be partially achieved is by comparing algebraic expressions to see if they are equivalent, i.e. to see if they are derivable from a common template or function definition. Herein a string matching and replacement algorithm is used to compare two functions. Initially, each expression is converted to a “reduced” form by reference to a set of standard rules so as, for example, to ensure that terms comprising the same variables are combined. The variables names in each term are then replaced with symbols to produce, for each expression, a corresponding character string which, arranged in a predetermined order, are compared to determine equivalence of the expressions. To achieve this, a matrix M is formed with entries corresponding to the symbols and the matrix elements are incremented according to pre-set rules. The matrix is then processed to identify equivalence.

BACKGROUND OF THE INVENTION

[0001] The present invention relates to the determination of the equivalence of two algebraic functions, such as those used in computer programs. The present invention includes a method and apparatus for carrying out such a determination and also to a computer program product including a computer readable medium having recorded thereon a computer program for performing such determination.

[0002] The present invention has, for its object, a technique which is specifically relevant to the determination of whether or not two algebraic expressions share a common template or function definition. Expressed mathematically, the objective is applicable as follows.

[0003] Given an algebraic expression f(x₁, x₂, . . . , x_(n)) of n-variables denoted by the symbol set {x₁, x₂, . . . , x_(n)} and another expression g(y₁, y₂, . . . , y_(n))of n-variables denoted by the symbol set {y₁, y₂, . . . , y_(n)} determine if there exists a permutation σ of degree n on the set {y₁, y₂, . . . , y_(n)} such that

g(y ₁ , y ₂ , . . . , y _(n))≡f(σ(y ₁), σ(y ₂), . . . , σ(y _(n))), n>0

where x ₁=σ(y ₁), x ₂=σ(y ₂), . . . , x _(n)=σ(y _(n)).

[0004] An identity would imply that f( ) and g( ) share a common template or function definition as mentioned earlier. The term “permutation” is defined as follows.

[0005] Let Ω be a set consisting of n distinct objects, which we may denote either by digits or by letters, say

Ω={1,2, . . . , n}.

[0006] Then a permutation on Ω implies the mapping of objects 1 to n to the same objects but in a different order—this is called mapping of Ω onto itself. A permutation σ of degree n is an operation on the set Ω which maps the digit i to the digit σ(i) in such a way that σ(i)≠σ(j) when i≠j. We write $\sigma = \begin{pmatrix} 1 & 2 & \ldots & \ldots & n \\ {\sigma (1)} & {\sigma (2)} & \ldots & \ldots & {\sigma (n)} \end{pmatrix}$

[0007] where the numbers σ(1), σ(2), . . . , σ(n) are an arrangement of the numbers 1, 2, . . . , n.

[0008] Note that since there are n! distinct arrangements possible of n objects, there are n! distinct permutations in the form given above. For example, the permutation of degree 5 given by $\sigma = \begin{pmatrix} 1 & 2 & 3 & 4 & 5 \\ 3 & 5 & 4 & 1 & 2 \end{pmatrix}$

[0009] in which σ(1)=3, σ(2)=5 . . . , σ(5)=2, may equally well be written in the form $\sigma = \begin{pmatrix} 4 & 2 & 3 & 1 & 5 \\ 1 & 5 & 4 & 3 & 2 \end{pmatrix}$

SUMMARY OF THE INVENTION

[0010] The invention comprises a method, apparatus and program product for determining the equivalence of two algebraic functions. In the invention, an algorithm is involved which comprises determining the symbol sets of the two functions; forming a matrix having a row for each symbol of one function and a column for each symbol of the other function; subjecting the functions to predetermined tests and incrementing the elements in the matrix in accordance with the results of said tests; and testing the matrix against predetermined rules for the symmetry of the matrix.

[0011] The algorithm may include comparing said symbol sets to determine whether they contain the same numbers of symbols.

[0012] The invention further provides a method, program and apparatus for determining the equivalence of two algebraic functions f( ) and g( ) having respective symbol sets {x₁, x₂, . . . , x_(n)} and {y₁, y₂, . . . , y_(n)} having equal numbers of symbols, the method comprising the steps of forming a matrix M with n rows and n columns and made up of elements M_(ij) corresponding to respective pairs of said symbols x_(i) and y_(j) where M_(ij) is calculated by initialising the values of M_(ij) to zero; putting x_(i)=0 in f( ) and y_(j)=0 in g( ) and replacing each other variable in f( ) and g( ) with a predetermined string character, then arranging the terms of f( ) and g( ) in ASCII order and comparing them and, if they match, incrementing the current value of M_(ij); putting x_(i)=1 in f( ) and y_(j)=1 in g( ) and replacing each other variable in f( ) and g( ) with a predetermined string character, then arranging the terms of f( ) and g( ) in ASCII order and comparing them and, if they match, incrementing the current value of M_(ij); putting x_(i)=−1 in f( ) and y_(j)=−1 in g( ) and replacing each other variable in f( ) and g( ) with a predetermined string character, then arranging the terms of f( ) and g( ) in ASCII order and comparing them and, if they match, incrementing the current value of M_(ij); ensuring that every element M_(ij) is either a zero or a three and that each row of the matrix M contains at least one three; for each entry M_(ij) which has the value three, exchanging the contents of columns i and j one with another; and inspecting the matrix M with the exchanged columns to determine if it is symmetrical.

[0013] The method and apparatus may further comprise analysing computer code containing first and second algebraic functions f(x₁, x₂, . . . , x_(n)) and g(y₁, y₂, . . . , y_(n)), the method being operable for analysing a relationship between said first and second functions and comprising the steps of:

[0014] (a) in each term in each function, replacing each variable name in that term with a predetermined symbol to form two character strings F( ) and G( ) each composed of sub-strings corresponding to respective terms of the respective function;

[0015] (b) arranging the sub-strings of each string in a like predetermined order; and

[0016] (c) comparing the strings one with another to determine equivalence of the two functions.

[0017] Preferably, the method further includes taking each variable x_(i) of the first function f(x₁, x₂, . . . , x_(n)) and running through variables y_(j) of the second function g(y₁, y₂, . . . , y_(n)) one-by-one and carrying out the steps of:

[0018] (a) removing all the terms in the first function which contain x_(i) and all the terms in the second function which contain y_(j) and, in each remaining term of each function, replacing each variable with a predetermined symbol and rearranging said terms in a predetermined order to form two character strings f₀( ) and g₀( );

[0019] (b) replacing each variable x_(i) and y_(j) by the value 1, replacing each remaining variable with a single predetermined symbol, and rearranging the terms of the functions in a predetermined order as two character strings f₁( ) and g₁( );

[0020] (c) replacing each variable x_(i) and y_(j) by the value −1, replacing each remaining variable with a single predetermined symbol, and rearranging the terms of the functions in a predetermined order as two character strings f⁻¹( ) and g⁻¹( );

[0021] (d) carrying out respective character string matching operations to determine if f₀( ) is identical to g₀( ), if f₁( ) is identical to g₁( ), and if f⁻¹( ) is identical to g⁻¹( ).

BRIEF DESCRIPTION OF THE DRAWINGS

[0022] For a better understanding of the invention, reference will be made, by way of example, to the accompanying drawings, in which:

[0023]FIG. 1 is a simplified diagram of a computer system; and

[0024]FIGS. 2A, 2B and 3 are flow-charts for explaining the method of the invention.

DETAILED DESCRIPTION

[0025] The invention comprises a method for determining the equivalence of two algebraic functions. The method may be implemented as a program for a computer, and the program may be stored on a storage medium, for example a CD-ROM, to form a program product according to the invention. Alternatively, a program product according to the invention may comprise a program made available for downloading from another computer. The computer program can be loaded into or made available to a suitable computer to form a computer system of the invention. FIG. 1 shows one embodiment of such a computer system.

[0026] This embodiment comprises a so-called stand-alone computer, i.e. one that is not permanently linked to a network. It includes a display monitor 2, a keyboard 3, a microprocessor—based central processing unit 4, a hard-disc drive 5 and a random access memory 6 all coupled one to another by a connection bus. The keyboard 3 is operable for enabling the user to enter commands into the computer along with user data. As well as keyboard 3, the computer may comprise, a mouse or tracker ball (not shown) for entering user commands especially if the computer is controlled by an operating system with a graphical user interface.

[0027] To load program instructions into the memory 6 and/or store them on the disc drive 5 so that the computer begins to operate or to become operable in accordance with the present invention, the computer 1 comprises a CD-ROM drive 8 for receiving a CD-ROM 9.

[0028] The program instructions are stored on the CD-ROM 9 from which they are read by the drive 8. However, as will be well understood by those skilled in the art, the instructions as read by the drive 8 may not be usable directly from the CD-ROM 9. Instead, they may be loaded into the memory 6 and stored in the hard disc drive 5 and used by the computer 1 from there. Also, the instructions may need to be decompressed from the CD-ROM using appropriate decompression software on the CD-ROM or in the memory 6 and may, in any case, be received and stored by the computer 1 in a sequence different to that in which they are stored on the CD-ROM.

[0029] In addition to the CD-ROM drive 8, or instead of it any other suitable input means could be provided, for, example a floppy-disc drive or a tape drive or a wireless communication device, such as an infrared receiver (none of these devices being shown).

[0030] The computer 1 also comprises a telephone modem 10 through which the computer is able temporarily to link up to the Internet via telephone line 11, a modem 12 located at the premises of an Internet service provider (ISP), and the ISP's computer 13. Also connected to the Internet are many remote computers, such as the computer 14, from which information, software and other services are available for downloading to the computer 1. Furthermore, the computer 1 does not have to be in a stand-alone environment. Instead, it could form part of a network (not shown) along with other computers to which it is connected on a permanent basis. It could also be permanently coupled to or have a temporary link to an intranet. An intranet is a group of data holding sites similar to Internet sites and arranged in the same way as the Internet but accessible only to particular users, for example the employees of a particular company. Instead of modem 10, the computer 1 could have a digital hard-wired link to the ISP's computer 13 or the computer 1 could itself comprise a permanently connected Internet site whether or not acting as an ISP for other remote users. Instead of the invention being usable only through the local keyboard 3, it may be available to remote users working through temporary or permanent link to computer 1 acting as ISP or simply as an Internet site.

[0031] Thus, instead of being provided by a local device such as the CD-ROM drive 8, the program instructions could be received or made available via the Internet from a remote computer such as the computer 14. Alternatively, the instructions could be received or made available from a computer connected with computer 1 via a network such as an intranet.

[0032] One embodiment of the method of this invention comprises a template matching algorithm to decide if two given algebraic expressions are derivable from a common expression template. Its possible applications include compiler-initiated optimisation of programs with common expression templates, and in the development of program/theorem proving techniques.

[0033] Thus, to perform the algorithm, the computer 1 is loaded with a suitable operating system, a compiler including a section for carrying out the equivalence determining algorithm, or simply a specific utility for carrying out that algorithm, and a section of source code which, it is assumed, contains the two expressions of which the equivalence is to be determined. The compiler/utility is supplied to computer 1 from a portable storage medium such as a series of floppy discs (not shown) or a CD-ROM or from another computer such as the computer 14. The algorithm will now be described by reference to the two flow charts of FIGS. 2A, 2B (FIG. 2B being a continuation of FIG. 2A) and 3. The numbered “steps” noted herein correspond to the step numbers of the flow charts.

[0034] Referring first to FIG. 2A, the two expressions f( ) and g( ) to be compared are identified in the relevant source code and then:

[0035] Step 1: Determine the Ordered Variables Symbol Sets

[0036] The expressions are parsed to identify the respective symbol sets and, in each set, the symbols are arranged in ASCII order (either ascending or descending). Say that the symbol set for f( ) is {x₁, x₂, . . . , x_(n)} and the set for g( ) is {y₁, y₂, . . . , y_(m)}.

[0037] Step 2: Check for Same Number of Symbols in Each Set

[0038] The respective numbers of symbols n and m in the two sets are compared and, if they are not the same, i.e. if n is not equal to m, the algorithm is exited, i.e. the method is abandoned. This is because the two expressions obviously cannot then be equivalent, i.e. they cannot then represent the same function. If n and m are the same, the method proceeds to Step 3.

[0039] Step 3: Reduce the Expressions

[0040] To facilitate the comparison of the two expressions, each is now reduced to a standardised format according to a predetermined set of rules.

[0041] This reduction operation, i.e. Step 3 of the algorithm, comprises a series of sub-operations or steps which will now be described with reference to FIG. 3.

[0042] Step 3.1: Make Variable Names Compliant

[0043] Initially, it is assumed that the expression to be converted is syntactically correct and does not contain any blanks. In the preferred embodiment, variable names appearing in the expression may comprise only lower-case alphabet characters, the underscore character, and digits. Also, a variable name may not start with a digit or end with an underscore. If these construction rules are not met, then the affected variable names are mapped (aliased) to alternative, but distinct, names obeying the construction rules, and these new names used instead.

[0044] Step 3.2: Remove Brackets

[0045] Brackets, if present, in the expression must be removed by carrying out the necessary operations needed to remove them, such as multiplying two parenthesised factors, discarding superfluous brackets, and so on.

[0046] Step 3.3: Insert Leading Operator if None Present

[0047] The expression is put in the following form:

[0048] <unitary operator><operand><operator><operand> . . . . <operator><operand>

[0049] where the unitary operator is either + (plus) or − (minus), and each operator is one of + (plus), − (minus), * (multiplication) or / (division). In the event that an expression does not commence with a unitary operator, one is added, i.e. a unitary operator + (plus) is inserted at the start of the expression. For example:

a*b/c becomes +a*b/c

[0050] Step 3.4: Map Division by Variable to Multiplication by Reciprocal of Variable

[0051] Division by a variable, for example the operator-operand pair /x, is replaced by multiplication by the reciprocal of the variable, where the reciprocal of the variable is formed as a new variable by appending an underscore to the variable, i.e./x is replaced by *x_ in the case of the given example.

[0052] Step 3.5: Map Division by Constant to Multiplication by Reciprocal of Constant

[0053] Division by a constant, for example /5, is replaced by multiplication by the reciprocal of the constant, i.e. /5 is replaced by *0.2 in the case of the given example.

[0054] Step 3.6: Replace “+” and “−” by Strings “+1*” and “−1*”

[0055] Next all + (plus) operators are substituted with the string “+1*” so that “+” becomes “+1*”. Similarly, all − (minus) operators are substituted with the string “−1*” so that “−” becomes “−1*”. Thus, for example:

+x becomes +1*x

and

−x*y+z becomes −1*x*y+1*z

[0056] Step 3.7: Convert Constants to e-format

[0057] Next the operands, which are constants (including the 1's introduced in the previous step) are converted into an e-format as follows:

“.[unsigned number]e[e-sign][unsigned exponent]”

[0058] where: [unsigned number] is a n-digit number comprising only digits and n is a predetermined fixed integer greater than 0; [e-sign] is the sign of the exponent and is one of > for plus or < for minus; and [unsigned exponent] is a m-digit number comprising only digits and m is a predetermined fixed integer greater than 0.

[0059] Thus, for example:

25=0.25*10² becomes 0.250000e>02

and

0.025=0.25*10⁻¹ becomes 0.250000e<01

[0060] where we have assumed n=6 and m=2. It is noted that any constant will be represented by a string of constant length m+n+3 characters in the e-format. Here e[e-sign] [unsigned exponent] represents the quantity 10 raised to the power [e-sign] [unsigned exponent], which must be multiplied to the number represented by. [unsigned number] to get the actual constant.

[0061] Now, the expression is free of the division operator and will contain at least one operand which is a constant. Each term in this expression will therefore have the following form:

<unitary operator><operand><*><operand> . . . <*><operand>

[0062] where the unitary operator is either + (plus) or − (minus), and between two consecutive operands is the multiplication operator *. After the terms are identified, the [e-sign] of each constant is restored from < or > to − or + respectively.

[0063] Step 3.8: Sort Operands in Each Term

[0064] In each term the operands are sorted (rearranged) in ascending order according to their ASCII value. This rearrangement is entirely permissible because the multiplication operator is commutative, i.e. the exchange of operands does not affect the result.

[0065] It is noted that no other variable will be able to place itself in the rearrangement between any particular variable and its reciprocal if they are both present. For example, if the variable “a” and its reciprocal “a_” are both present, they will be sorted so as to remain together as “*a*a_”.

[0066] Step 3.9: Eliminate Variable/Reciprocal Pairs

[0067] Next, all operator-operand sequences of the form “*a*a_” are eliminated from the term. An expression, such as, a³/a² will appear as “*a*a*a*a_*a_”. After “*a*a_” has been eliminated from it, “*a*a*a_” will remain, from which “*a*a_” must, again, be eliminated. That is, the elimination process must be continued till no further elimination is possible.

[0068] Step 3.10: Consolidate Constants in Each Term

[0069] After the sorting step 3.8, the operands, which are constants in the term will be bunched up at the beginning of the terms where they can be easily identified and replaced by a single constant. Thus, for example:

+0.100000e+01*s*k*m*s_(—)*0.500000e+00

[0070] after arranging the operands in ascending order becomes

+0.100000e+01*0.500000e+00*k*m*s*s_(—)

[0071] and after eliminating the variables s and s_ and consolidating the constants, the term becomes

+0.500000e+00*k*m

[0072] At this stage a term will have the following form:

<unitary operator><constant><*><operand> . . . <*><operand>

[0073] where each operand is a variable name, whose ASCII value is not lower than that of its preceding operand, if any. This is the reduced form of a term. In the reduced form, the non-constant part of a term is called a variable-group. For example, if the term in the reduced form is “+0.250000e+01*m*m*s”, then its variable-group is “*m*m*s”.

[0074] Step 3.11: Consolidate Like Terms

[0075] In an expression, all those terms whose variable-groups match, are combined by modifying the constant in one of the terms, and eliminating the others.

[0076] Step 3.12: Sort Terms in the Expression

[0077] Finally, the reduced terms in the expression are rearranged in an ascending order according to the ASCII value of their respective variable-groups. In this final form, the expression is said to be in its reduced form. Note, in particular, that no two terms in a reduced expression will have the same variable-group.

[0078] Henceforth, references to the expressions f( ) and g( ) refer to the reduced form of those expressions.

[0079] Returning now to FIG. 2A and then its continuation FIG. 2B, the analysis algorithm then continues with the following steps:

[0080] Step 4: Are Variables Symbols Sets Identical?

[0081] If the variables in the two ordered sets {x₁, x₂, . . . , x_(n)} and {y₁, y₂, . . . , y_(n)} are identical (that is, for each i=1 to n, x_(i) and y₁ have identical symbols) then a simple string match between the reduced strings of f( ) and g( ) will tell us if f( )=g( ) or not. If the strings match then f( )=g( ), else not; and the algorithm ends. If the variable sets are not identical, g( ) to Step 5.

[0082] Step 5: Is n=1?

[0083] If n is not equal to 1 go to Step 8.

[0084] Step 6: Put y₁=x₁

[0085] Put y₁=x₁ in g( ) to form g(x₁).

[0086] Step 7: Compare f(x₁) and g(x₁)

[0087] If f(x₁)=g(x₁) exit with condition true, else exit with condition false. Steps 6 and 7 only apply to the case where n=1, i.e. where two expression each have only one variable. The comparison indicates whether or not the expressions are the same and the method then terminates.

[0088] Step 8: Reference Variables

[0089] This and the following steps are carried out where n is not equal to 1. Replace each of the variables in f( ) with the symbol #, and rearrange the terms (sign included), treating them as character strings, in ascending (or, descending) order. Call this new string F( ). Likewise create G( ) from g( ).

[0090] Step 9: Compare F( ) and G( )

[0091] If the strings F( ) and G( ) do not match, then f( )≠g( ) and the algorithm ends, i.e. exit with condition false. If they match, g( ) to the next step. Note that replacing each variable by the symbol # makes all variables have the same character count. Hence if the strings F( ) and G( ) do not match, it can be definitely concluded that f( ) and g( ) cannot be identical. However, if they match, we cannot yet conclude that f( ) and g( ) are identical.

[0092] Step 10: Set Up Matrix M

[0093] Set up a square table of n rows and n columns, with the columns labelled y₁, y₂, . . . , y_(n), and the rows labelled x₁, x₂, . . . , x_(n). Call this matrix M. Let an element in this matrix corresponding to x_(i) and y_(j) be designated by M_(ij). Initialise all the elements in the matrix M to 0 (zero).

[0094] Step 11: Calculate Elements of Matrix M

[0095] To calculate the value of M_(ij) do the following three calculations:

[0096] (a) Put x_(i)=y_(j)=0 in f( ) and g( ), that is, remove all the terms in f( ) which contain x_(i) and all the terms in g( ) which contain y_(j). Replace the remaining variables in f( ) and g( ) with the symbol #, and rearrange their respective terms (sign included), treating them as character strings, in ascending (or, descending) order. Call the resulting strings of f( ) and g( ), respectively, f₀( ) and g₀( ). If the strings f₀( ) and g₀( ) match, increment the current value of M_(ij) by 1.

[0097] (b) Put x_(i)=y_(j)=1 in f( ) and g( ), that is, remove all occurrences of x_(i) and its immediately preceding multiplication operator from f( ) and do likewise for all occurrences of y_(j) in g( ). Replace the remaining variables in f( ) and g( ) with the symbol #, and rearrange their respective terms (sign included), treating them as character strings, in ascending (or, descending) order. Call the resulting strings of f( ) and g( ), respectively, f₁( ) and g₁( ). If the strings f₁( ) and g₁( ) match, increment the current value of M_(ij) by 1.

[0098] (c) Put x_(i)=y_(j)=−1 in f( ) and g( ), that is, remove x_(i) and its immediately preceding multiplication operator from all the terms in f( ) which contain x_(i) for each such term change its sign if an odd number of x_(i)s were removed from it, and do likewise for all occurrences of y_(j) in g( ). Replace the remaining variables in f( ) and g( ) with the symbol #, and rearrange their respective terms (sign included), treating them as character strings, in ascending (or, descending) order. Call the resulting strings of f( ) and g( ), respectively, f⁻¹( ) and g⁻¹( ). If the strings f⁻¹( ) and g⁻¹( ) match, increment the current value of M_(ij) by 1.

[0099] Step 12: If Any Element of M is Not 0 or 3, Exit

[0100] If the matrix M has elements, whose values are other than 0 or 3, then f( )≠g( ), and the algorithm ends, i.e. exit with condition false. This is because if f( ) and g( ) are equal then all the tests, a, b and c of Step 11, together should either fail or succeed for a given (x_(i), y_(j)) pair. If the elements of M are either 0 or 3, go to the next step.

[0101] As an explanation of step 12, if two functions f( ) and g( ) are equivalent such that the variable x_(i) in f( ) can correspond to the variable y_(j) in g( ), then their respective strings f₀( ) and g₀( ) derived from test a of step 11 must obviously be identical to each other. Likewise for tests b the strings f₁( ) and g₁( ) must be identical and for test c the strings f⁻¹( ) and g⁻¹( ) must be identical. Thus, the element of M corresponding to these tests must be 1+1+1=3. Similarly, if f( ) and g( ) are equivalent, but x_(i) in f( ) cannot correspond to y_(j) in g( ), then all tests a, b, c in step 11 will fail and the element of M corresponding to these tests will be 0. It is important to note that these tests are carried out on f( ) and g( ) only when each of them is in its reduced form prior to the tests and F( )=G( ). This fact is important and forms the basis for making the string comparisons.

[0102] When two functions are not equivalent they can never lead to a matrix element whose value is 3 because at least one of the steps a, b, or c will result in a string mismatch. They may, however, lead to a matrix element whose value is 0. See Examples 2 and 3.

[0103] Below, 3 examples are provided to illustrate the calculation of the elements of M to support the above noted point.

EXAMPLE 1

Let f(x, y)≡x+x*y and g(a, b)≡a+a*b

[0104] where it is clear that if x=a and y=b, then the functions f and g are equivalent, but not if x=b and y=a.

[0105] Consider the case when x=a when the functions are equivalent.

[0106] Following step 11(a) we put x=a=0 and get

f ₀( )≡0; g ₀( )≡0.

[0107] Clearly, the strings of f₀( ) and g₀( ) will match and M_(ij)=1, where i refers to variable x and j refers to variable a.

[0108] Following step 11(b) we put x=a=1 and get

f ₁( )≡1+y; g ₁( )≡1+b.

[0109] After substituting y and b with the character ‘#’ and completing the remaining substeps of step 11(b) we will once again find that the strings of f₁( ) and g₁( ) will match and M_(ij)=2.

[0110] Following step 11(c) we put x=a=−1 and get

f ⁻¹( )≡−1−y; g ⁻¹( )≡−1−b.

[0111] After substituting y and b with the character ‘#’ and completing the remaining substeps of step 11(c) we will once again find that the strings of f⁻¹( ) and g⁻¹( ) will match and M_(ij)=3.

[0112] Similarly, the element of M corresponding to i=y and j=b will also evaluate to 3.

[0113] On the contrary, the element of M corresponding to x=b will evaluate to 0, as follows:

[0114] Following step 11(a) we put x=b=0 and get

f ₀( )≡0; g ₀( )≡a.

[0115] Clearly, the strings of f₀( ) and g₀( ) will not match and M_(ij)=0, where i refers to variable x and j refers to variable b.

[0116] Following step 11(b) we put x=b=1 and get

f ₁( )≡1+y; g ₁( )≡a+a.

[0117] After completing the remaining substeps of step (b) we will once again find that the strings of f₁( ) and g₁( ) will not match and M_(ij)=0.

[0118] Following step 11(c) we put x=b=−1 and get

f ⁻¹( )≡−1−y; g ⁻¹( )≡a−a.

[0119] After completing the remaining substeps of step (c) we will once again find that the strings of f⁻¹( ) and g⁻¹( ) will not match and M_(ij)=0.

[0120] Similarly, the element of M corresponding to i=y and j=a will also evaluate to 0.

EXAMPLE 2

Let f(x, y)≡x+x*y+z and g(a, b)≡a+b*c+a

[0121] where it is clear that the functions f( ) and g( ) are not equivalent. When two functions are not equivalent they can never lead to a matrix element whose value is 3.

[0122] Consider the case when x=a.

[0123] Following step 11(a) we put x=a=0 and get

f ₀( )≡z; g ₀( )≡b*c.

[0124] Clearly, the strings of f₀( ) and g₀( ) will not match and M_(ij)=0, where i refers to variable x and j refers to variable a.

[0125] Following step 11(b) we put x=a=1 and get

f ₁( )≡1+y+z; g ₁( )≡1+b*c+1.

[0126] After completing the remaining substeps of step 11(b) we will find that the strings of f₁( ) and g₁( ) will not match and M_(ij)=0.

[0127] Following step 11(c) we put x=a=−1 and get

f ⁻¹( )≡−1−y+z; g ⁻¹( )≡−1+b*c−1.

[0128] After completing the remaining substeps of step 11(c) we will find that the strings of f⁻¹( ) and g⁻¹( ) will not match and M_(ij)=0, which is also the final value for this element.

[0129] Now consider the case when x=b. Following step 11(a) we put x=b=0 and get

f ₀( )≡z; g ₀( )≡a+a.

[0130] Clearly, the strings of f₀( ) and g₀( ) will not match and M_(ij)=0, where i refers to variable x and j refers to variable b.

[0131] Following step 11(b) we put x=b=1 and get

f ₁( )≡1+y+z; g ₁( )≡a+c+a.

[0132] After completing the remaining substeps of step 11(b) we will find that the strings of f₁( ) and g₁( ) will not match and M_(ij)=0.

[0133] Following step 11(c) we put x=b=−1 and get

f ⁻¹( )≡−1−y+z; g ⁻¹( )≡a−c+a.

[0134] After completing the remaining substeps of step 11(c) we will find that the strings of f⁻¹( ) and g⁻¹( ) will not match and M_(ij)=0, which is also the final value for this element.

[0135] Similarly, other elements of M will all show that none of them have the value 3.

EXAMPLE 3

Let f(x, y)≡x+x*y+z and g(a, b)≡a+b+c

[0136] where it is clear that the functions f( ) and g( ) are not equivalent. When two functions are not equivalent they can never lead to a matrix element whose value is 3.

[0137] Consider the case when x=b.

[0138] Following step 11(a) we put x=b=0 and get

f ₀( )≡z; g ₀( )≡a+c.

[0139] Clearly, the strings of f₀( ) and g₀( ) will not match and M_(ij)=0, where i refers to variable x and j refers to variable b.

[0140] Following step 11(b) we put x=b=1 and get

f ₁( )≡1+y+z; g ₁( )≡a+1+c.

[0141] After completing the remaining substeps of step 11(b) we will find that the strings of f₁( ) and g₁( ) will match and M_(ij)=1.

[0142] Following step 11(c) we put x=b=−1 and get

f ⁻¹( )≡−1−y+z; g ⁻¹( )≡a−1+c.

[0143] After completing the remaining substeps of step 11(c) we will find that the strings of f⁻¹( ) and g⁻¹( ) will not match and M_(ij)=1, which is also the final value for this element.

[0144] Similarly, other elements of M will all show that none of them have the value 3 but they need not be calculated since even a single element whose value is neither 0 nor 3 in M is sufficient to show that the two functions are not equivalent.

[0145] The three steps (a), (b), (c) are heuristics. They were chosen because step (a) will remove terms containing a given variable, step (b) will delete a given variable from terms containing it, step (c), in addition to what step (b) does, will change the sign of terms which contain a given variable in odd multiples. The reduced strings of f( ) and g( ) will be affected differently in at least in one of the steps if f( ) is not equal to g( ), assuming that F( ) and G( ) were equal. These three steps appear to be sufficient for deciding whether f( ) and g( ) are equal or not. However, as a precaution, finally, f( ) and g( ) are compared with the appropriate substitution of variables given by a permutation to ensure that the result is indeed correct.

[0146] Step 13: If Any Row of M Does Not Contain 3, Exit

[0147] For each row, beginning with the i=1 row, scan the columns of the matrix for a 3. If the row is devoid of a 3, then f( )≠g( ), and the algorithm ends, i.e. exit with condition false.

[0148] Step 14: Exchange Columns of M

[0149] For each row, beginning with i=1 row, check the row's columns. For each column j>i in the ith row for which M_(ij)=3 exchange the contents of columns i and j in M (including the column labels). Note that the column numbering is not changed, only the contents of the columns and their labels are exchanged. Note also that the exchange operations amount to finding the permutation σ of degree n on the set {y₁, y₂, . . . , y_(n)}.

[0150] Step 15: If M is Not Symmetric, Exit

[0151] Examine the resulting matrix M. If M is not symmetric then f( )≠g( ), and the algorithm ends, i.e. exits with condition false. If it is symmetric then go to Step 16.

[0152] Step 16: Set Up the Permutation σ

[0153] Set up the permutation σ, that is a correspondence between the elements of the variables set {x₁, x₂, . . . , x_(n)} and the variables set {y₁, y₂, . . . , y_(n)} on the basis of the diagonal of M. This can be done by going down the diagonal of M and, for each diagonal element, creating the variable labels pair (x_(i), y_(k)), where y_(k) is the current column label of the i-th column.

[0154] Step 17: Compare f( ) and g( )

[0155] In g( ), now replace each occurrence of the variable y_(k) by its corresponding x_(i) for k=1 to n or, alternatively replace each occurrence of the variable x_(i) by y_(k) for i=1 to n in f( ). If in their reduced form, f( ) and g( ) are equivalent, then exit from the algorithm with condition true, else exit with condition false.

[0156] Step 18: Off-diagonal Elements with Value 3

[0157] If M has off-diagonal elements, which are 3s, it indicates that alternative permutations σ of degree n on the set {y₁, y₂, . . . , y_(n)} exist. This situation arises if f( ), and therefore g( ), is symmetric with respect to two or more of its variables. Let us call a column a 3-column if it has an off-diagonal 3 in any of its rows. An alternative permutation is obtained by exchanging the contents of such a 3-column with its symmetric counterpart on the other side of the diagonal (that is, if the 3-column chosen is the j-th column then its symmetric counterpart will be the (n−j)-th column).

FURTHER EXAMPLES EXAMPLE 4

[0158] Let

f(x, y, z)≡+0.100000e+01*x+0.400000e+01*x*y+0.200000e+01*x*y*z+0.300000e+01*y*z

and

g(a, b, c)≡+0.100000e+01*b+0.400000e+01*b*c+0.200000e+01*a*b*c+0.300000e+01*a*c

[0159] (The blank spaces on the right hand side of the identities have been introduced for visual clarity here, the reduced form does not have them.) By inspection we can deduce that if x=b, y=c, z=a, then the two functions are indeed identical. Let us see if the algorithm described above will also provide the same conclusion.

[0160] Step 4.

[0161] Since the symbol sets {x, y, z} and {a, b, c} are different, move to Step 5.

[0162] Step 5.

[0163] Since n≠1, move to Step 8.

[0164] Step 8.

[0165] Substitute every variable in f(x, y, z) and g(a, b, c) with the character # to obtain

F≡+0.100000e+01*#+0.400000e+01*#*#+0.200000e+01*#*#*#+0.300000e+01*#*#

and

G≡+0.100000e+01*#+0.400000e+01*#*#+0.200000e+01*#*#*#+0.300000e+01*#*#

[0166] Rearrange the terms in ascending order to obtain

F≡+0.100000e+01*#+0.200000e+01*#*#*#+0.300000e+01*#*#+0.400000e+01*#*#

and

G≡+0.100000e+01*#+0.200000e+01*#*#*#+0.300000e+01*#*#+0.400000e+*#*#

[0167] Step 9.

[0168] A string match between F and G will show that they match, so move to next step.

[0169] Step 10.

[0170] The initial matrix M a b c x 0 0 0 y 0 0 0 z 0 0 0

[0171] Step 11.

[0172] The final matrix M a b c x 0 3 0 y 0 0 3 z 3 0 0

[0173] Step 12.

[0174] The matrix M has no elements which are neither 0 or 3, so move to Step 13.

[0175] Step 13.

[0176] Each row of matrix M has a 3 so go to Step 14.

[0177] Step 14.

[0178] Exchange columns.

[0179] First and second columns b a c x 3 0 0 y 0 0 3 z 0 3 0

[0180] Second and third columns b c a x 3 0 0 y 0 3 0 z 0 0 3

[0181] Steps 15 to 17.

[0182] Note that the matrix M is diagonal (and hence also symmetric), and the diagonal relates the row labels with the column labels as x=b, y=c, z=a. If we substitute for a, b, c, respectively, by z, x, y in g( ), and use the reduction algorithm described in connection with FIG. 3, we will find that f( )≡g( ).

EXAMPLE 5

[0183] Let

f(x, y, z)≡+0.100000e+01*x*x+0.100000e+01*x*y+0.100000e+01*y*y+0.100000e+01*z

and

g(a, b, c)≡+0.100000e+01*a*a+0.100000e+01*a*b+0.100000e+01*b*b+0.100000e01*c

[0184] By inspection we can deduce that if x=a, y=b, z=c, then the two functions are indeed identical. Alternatively, if x=b, y=a, z=c, then again the two functions are identical. Let us see if the algorithm described above will provide the same conclusions.

[0185] Step 4.

[0186] Since the symbol sets {x, y, z} and {a, b, c} are different, move to Step 5.

[0187] Step 5.

[0188] Since n≠1, move to Step 8.

[0189] Step 8.

[0190] Substitute every variable in f(x, y, z) and g(a, b, c) with the character # to obtain

F≡+0.100000e+01*#*#+0.100000e+01*#*#+0.100000e+01*#*#+0.100000e01*#

and

G≡+0.100000e+01*#*#+0.100000e+01*#*#+0.100000e+01*#*#+0.100000e01*#

[0191] Rearrange the terms in ascending order to obtain

F≡+0.100000e+01*#+0.100000e+01*#*#+0.100000e+01*#*#+0.100000e+01*#*#

and

G≡+0.100000e+01*#+0.100000e+01*#*#+0.100000e+01*#*#+0.100000e+01*#*#

[0192] Step 9.

[0193] A string match between F and G will show that they match, so move to Step 10.

[0194] Step 10.

[0195] The initial matrix M a b c x 0 0 0 y 0 0 0 z 0 0 0

[0196] Step 11.

[0197] The final matrix M a b c x 3 3 0 y 3 3 0 z 0 0 3

[0198] Step 12.

[0199] The matrix M has no elements which are neither 0 or 3, so move to Step 13.

[0200] Step 13.

[0201] Each row of matrix M has a 3 so go to Step 14.

[0202] Step 14.

[0203] Since the matrix M already appears in symmetric form with all its diagonal elements containing a 3, exchange of columns is implicitly complete.

[0204] Steps 15 to 18.

[0205] Note that in matrix M the diagonal relates the row labels with the column labels as x=a, y=b, z=c. If we substitute for a, b, c, respectively, by x, y, z in g( ), again using the reduction algorithm described with reference to FIG. 3 of the drawings, it will be found that f( ) g( ). Also notice that we can exchange columns 1 and 2 and obtain b a c x 3 3 0 y 3 3 0 z 0 0 3

[0206] And thence, x=b, y=a, z=c. So both the solutions are found by the algorithm.

[0207] It will be understood that the invention is not limited to the specific embodiments described above but includes modifications and developments within the purview of those skilled in the art and limited solely by the following claims. 

1. A method for determining the equivalence of two algebraic functions, the method including the steps of:— (a) determining the symbol sets of the two functions; (b) forming a matrix having a row for each symbol of one function and a column for each symbol of the other function; (c) subjecting the functions to predetermined tests and incrementing the entries in the matrix in accordance with the results of said tests; and (d) testing the matrix against predetermined rules for the symmetry of the matrix.
 2. A method according to claim 1, including, between steps (a) and (b), carrying out the step of:— (e) comparing said symbol sets to determine whether they contain the same numbers of symbols.
 3. A method for determining the equivalence of two algebraic functions f( ) and g( ) having respective symbol sets {x₁, x₂, . . . , x_(n)} and {y₁, y₂, . . . , y_(n)} having equal numbers of symbols, the method comprising the steps of:— (a) forming a matrix M with n rows and n columns and made up of elements M_(ij) corresponding to respective pairs of said symbols x_(i) and y_(j) where M_(ij) is calculated by: (i) initialising the values of M_(ij) to zero; (ii) putting x_(i)=0 in f( ) and y_(j)=0 in g( ) and replacing each other variable in f( ) and g( ) with a predetermined string character, then arranging the terms of f( ) and g( ) in ASCII order and comparing them and, if they match, incrementing the current value of M_(ij); (iii) putting x_(i)=1 in f( ) and y_(j)=1 in g( ) and replacing each other variable in f( ) and g( ) with a predetermined string character, then arranging the terms of f( ) and g( ) in ASCII order and comparing them and, if they match, incrementing the current value of M_(ij); (iv) putting x_(i)=−1 in f( ) and y_(j)=−1 in g( ) and replacing each other variable in f( ) and g( ) with a predetermined string character, then arranging the terms of f( ) and g( ) in ASCII order and comparing them and, if they match, incrementing the current value of M_(ij); (b) inspecting the matrix M to determine if every element M_(ij) is either a zero or a three and that each row of the matrix M contains at least one three; (c) for each element M_(ij) which has the value three and satisfies the relationship j>i, exchange the contents of columns i and j one with another; (d) inspecting the matrix M with the exchanged columns to determine if it is symmetrical.
 4. A method according to claim 1, the method being operable for analysing a relationship between said first and second functions and comprising the steps of:— (f) in each term in each function, replace each variable name in that term with a predetermined symbol to form two character strings F( ) and G( ) each composed of sub-strings corresponding to respective terms of the respective function; (g) arranging the sub-strings of each string in a like predetermined order; and (h) comparing the strings one with another to determine equivalence of the two functions.
 5. A method according to claim 4, including taking each variable x_(i) of the first function f(x₁, x₂, . . . , x_(n)) and running through variables y_(j) of the second function g(y₁, y₂, . . . , y_(n)) one-by-one and carrying out the steps of:— (i) removing all the terms in the first function which contain x_(i) and all the terms in the second function which contain y_(j) and, in each remaining term of each function, replacing each variable with a predetermined symbol and rearranging said terms in a predetermined order to form two character strings f₀( ) and g₀( ); (j) replacing each variable x_(i) and y_(j) by the value 1, replacing each remaining variable with a single predetermined symbol, and rearranging the terms of the functions in a predetermined order as two character strings f₁( ) and g₁( ); (k) replacing each variable x_(i) and y_(j) by the value −1, replacing each remaining variable with a single predetermined symbol, and rearranging the terms of the functions in a predetermined order as two character strings f⁻¹( ) and g⁻¹( ); (l) carrying out respective character string matching operations to determine if f₀( ) is identical to g₀( ), if f₁( ) is identical to g₁( ), and if f⁻¹( ) is identical to g⁻¹( ).
 6. A method according to claim 1, comprising the steps of:— (m) inspecting the variables in the sets of variables {x₁, x₂, . . . , x_(n)} and {y₁, y₂, . . . , y_(n)} of the respective expressions f( ) and g( ) to determine if they are identical; (n) if the variable sets are not identical, and if n=1, putting y₁=x₁ in g( ) and determining if the strings of f( ) and g( ) now match; (p) replacing each of the variables in f( ) and g( ) with the symbol #, and rearranging the terms (sign included) in ascending (or, descending) order to produce new strings F( ) and G( ).
 7. Apparatus for determining the equivalence of two algebraic functions, the apparatus comprising:— (i) means for determining the symbol sets of the two functions; (ii) means for forming a matrix having a row for each symbol of one function and a column for each symbol of the other function; (iii) means for subjecting the functions to predetermined tests and incrementing the elements in the matrix in accordance with the results of said tests; and (iv) means for testing the matrix against predetermined rules for the symmetry of the matrix.
 8. Apparatus according to claim 7, the apparatus comprising a processor and data memory having stored therein a computer program to be optimised and also computer code for causing the processor to be operable to carry out the following actions: (a) identifying within said program first and second algebraic functions f(x₁, x₂, . . . , x_(n)) and g(y₁, y₂, . . . , y_(n)); (b) in each term in each function, replacing each variable name x₁, x₂, . . . , x_(n) and y₁, y₂, . . . , y_(n) in that term with a single predetermined symbol to form two character strings F( ) and G( ) each composed of sub-strings corresponding to respective terms of the respective function; (c) arranging the sub-strings of each string in a like predetermined order; and (d) comparing the strings one with another to determine equivalence of the two functions.
 9. A program product comprising a body of software code for determining the equivalence of two algebraic functions, the software code including:— (a) a code portion for determining the symbol sets of the two functions; (b) a code portion for forming a matrix having a row for each symbol of one function and a column for each symbol of the other function; (c) a code portion for subjecting the functions to predetermined tests and incrementing the elements in the matrix in accordance with the results of said tests; and (d) a code portion for testing the matrix against predetermined rules for the symmetry of the matrix.
 10. A product according to claim 9, including a code portion for comparing said symbol sets to determine whether they contain the same numbers of symbols.
 11. A compiler comprising software code which is operable for carrying out the method of claim
 1. 